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<N , Abstract 

The microscopic structure and movement of reaction fronts in reaction diffusion systems far 
from equilibrium are investigated. We show that some three-site interaction models exhibit exact 
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diffusive shock measures, i.e. domains of different densities connected by a sharp wall without 
correlations. In all cases fluctuating domains grow at the expense of ordered domains, the absence 
of growth is possible between ordered domains. It is shown that these models give rise to aspects 
not seen in nearest neighbor models, viz. double shocks and additional symmetries. A classification 
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of the systems by their symmetries is given and the link of domain wall motion and a free fermion 
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I. INTRODUCTION 

The emergence of patterns and fronts is a challenging problem in biology, chemistry and 
physics, for a review see Ref. |l|. In biology for example bacteria aggregate building u 
regions with a high density in coexistence with regions with a low density of organisms 
A typical example in chemistry is the movement of reaction fronts. In physics the movement 
of domain walls is directly related to the problem of coarsening, for example in magnetic 
systems. In Ref. |3| the phenomenon of hysteresis in a driven diffusive system is explained 
by the movement of shocks, i.e., a jump in the density profile. Various other phenomena in 
many particle systems can be attributed to the emergence of shocks, for example the first 
order transition in the phase diagram of the asymmetric exclusion process (ASEP) Q or 
phase coexistence in a driven diffusive system coupled to reaction kinetics p,ly]- Recently, 
shocks in quantum systems have also been discussed [7|, |8|. 

On a macroscopic level shock fronts are described by partial differential equations. The 
most prominent equations are the Fisher and the Burgers' equations. The Fisher equation 

a 

d d 2 

— p(x, t) = D—p(x, t) + ap(x, t) - bp 2 (x, t) (1) 

was originally proposed as a model for the propagation of a mutant gene. It shows traveling 
wave solutions and may be used for modeling systems without conservation of the order 
parameter. 

The inviscid Burgers' equation [lXJ] 

-p(x, t) = -a— [p{x, t) (1 - p(x, t))] (2) 

was proposed as a model for turbulent fluid motion. It as well shows shock solutions, but it 
may be used for modeling systems with particle conservation. 

In this paper we demonstrate for some models how these macroscopic shocks originate 
from the microscopic dynamics. It is known that some driven diffusive systems can be 
described by the Fisher or Burgers' equation in the hydrodynamic limit. This limit is 
achieved by scaling the lattice constant to zero while keeping the overall length of the system 
constant, the time has to be rescaled appropriately. One of these models is the ASEP which 
is in the hydrodynamic limit described by the Burgers' equation. Another model is the 
branching and coalescing random walk, which is in the hydrodynamic limit described by the 



Fisher equation |ll|. In these two models the microscopic dynamics can be described by 
exact shock measures. An exact shock measure is a state where two product measures with 
different densities are connected. The time evolution of this state is given by a diffusion 
elation with respect to the position of the density step Q Q Q. In particndar this 
implies that the microscopic structure of the system is known at all times. 

The physical properties of large classes of one-species reaction-diffusion models with 
nearest neighbor interactions have been widely studied 14, jjl 16L Il7|. There are only four 
known models which show shocks without correlations 
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the ASEP, the branching 



and coalescing random walk (BCRW), the asymmetric Kawasaki-Glauber process (AKGP) 
and the brick layer model. While the former three models are exclusion models where the 
number of particles on each lattice site is restricted to at most one, the particle number 
is not restricted in the latter one. Here, the investigation shall be extended to three-site 
interaction exclusion m ode.s Q Q and it will be shown that three-site interactions give 
rise to models with exact shock measures which show aspects not seen in nearest neighbor 
models, viz. double shocks and additional symmetries even though no free fermion condition 
is met. 

At this point we remark that the question of phase separation is directly linked to the 
movement of domain walls because coarsening is generic for this phenomenon. But although 
some of the models presented in this paper show growing domains we will argue that this 
mechanism cannot be used for constructing nonequilibrium models with two species (one 



type of particles and vacancies) showing phase separation in one dimension [22J, |2i 

We present a survey of one-species models with three-site interactions and open bound- 
aries whose time evolution is described by an exact diffusive shock measure, to be defined 
below. For the cases of shocks between two nonfluctuating phases (densities to 1) and two 
fluctuating phases (both densities are different from and 1) the list is complete. For the 
case of a shock from a nonfluctuating phase (density or 1) to a fluctuating phase (density 
between and 1) the variety of models is too large to give a complete survey - the number 
of free parameters rises from 12 to 56 when the interaction range is increased from two to 
three. But we present a classification of models with respect to their symmetries, where we 
considered models where at least two of the symmetries charge (C), parity (P) and time 
(T) are valid independently. 

We also address the question to what extent the description of the dynamics of the model 
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by the movement of shock fronts is sufficient. To this end the interactions of the domain 
walls are determined and possible parallels to free fermion systems are discussed. 

II. FORMALISM 

On each lattice site k {k = 1, . . . , L) of a one-dimensional lattice there may be at most 
one particle (A) or a vacancy (0). One can also consider these two-state systems as spin 
systems, a particle is represented by a down spin (||)) and a vacancy by an up spin (|T))- The 
stochastic dynamics of the models are defined by a master equation [25] which is convenient! 
expressed in the quantum-Hamiltonian formalism for spin-1/2 chains as described in Ref 
To each lattice configuration 77 we assign a basis vector \rj) which is given by the tensor 
product of the single-site states. The probability vector describing the system can then be 
written as 

\P(t)) = J2 P (v,t)\v) (3) 

v 

where P(rj, t) is the probability at time t to find the system in the state rj. The time evolution 
of the system is described by a master equation which can be written as 

j t \P{t)) = -H\P{t)) (4) 

where H is the stochastic generator of the process. Due to the analogy of equation (J3J) to the 
Schrodinger equation (in imaginary time) H is often called the Hamiltonian of the system. 
Conservation of probability requires that the sum of the entries of each column is zero, 

(s\H = 0, (5) 

where (s\ = (1, 1, ...) is the so-called summation vector. 

For finite interaction range and spatially homogeneous kinetics it is convenient to write 
the generator as 

# = -$>k-&i-&L-i, (6) 

k 

where the local Hamiltonians hk contain the rates of the elementary local transitions and b\, 
&l_i account for events at the left respectively right boundary. The operators h m include 
only operators acting on the sites m, m + 1 and m + 2; the operators b m include only 
operators acting on m and m + 1. They are formulated using the particle number operator 



(n k ), vacancy number operator {y k — 1 — n k ), the particle creation (s^) and annihilation 
(s^) operators. The lower index represents the lattice site on which the respective operator 
acts. The diagonal entries of the Hamiltonian have to be chosen such that conservation of 
probability is fulfilled which can be easily constructed by considering 

(s\s~ = (s\v k ; (s\s+ = {s\n k . (7) 

For example diffusion to the right (A0 — ► 0A) with rate D r is written as h k = 
A- { s t s k+i -n k v k+1 ). 

III. PRODUCT AND SHOCK MEASURES 

In what follows |-) 1 denotes a probability vector for a single site and an operator without 
index a single site operator. If no correlations are present the probability vectors are simply 
given by 

\P(t)) = \p) = \p)i®-"®\p)i (8) 

with 

ip)i=piDi +(i-p)it>i = rj- w 

In this case the calculation of the action of the stochastic generator H can be simplified by 
using the following identities: 

s- \p), = l —^n \p), ; s+ \p), = -^-v \p) x . (10) 

A product measure \p) is a stationary state of the system, if 

H\p)=0. (11) 

As a shock measure \p\, p2,k) we define the state which is a product measure with density 
P\ up to and including site k and beginning from site fc + 1 a product measure with density 

\p 1 ,p 2 ,k) = \p 1 )f®\p 2 f 1 L - k . (12) 

We are interested in those systems, for which the time evolution of the shock measure is 



given by a diffusion equation (see figure HJ) 

— |pi, p 2 , k) = -H \pi,p 2 , k) = Si \pi,p 2 , k - 2) + 5 2 \pi, p 2 , k-1) 

+ $3 \pi, P2, k + 1) + 5a \pi, Pi, k + 2) 
- h\PuPz,k) . (13) 

Starting from an initial density step the system will evolve into an exact diffusive shock 
measure defined by 

\P(t))=J2Pi(t)\pnP2J) (14) 

i 

which is a time dependent superposition of sharp shocks weighted with pi(t). This can be 
seen by solving Eq. (J13)) for an initial shock at position k: 

\p h p 2 , k,t) =^2 G ( 1 ' f \ k i °) IPi> Pa, 
I 

1 f 77 

G(l,t\k,0) = — dpexp(-e p t + i(k-l)p) ( 15 ) 

27T J-7T 

e p = —(Si exp(-2ip) + 5 2 exp(-ip) + 5 3 exp(ip) + S 4 exp(2ip) - S 5 ). 

G(l,t\k,0) is the Green's function for this problem, for large arguments (k — I) and late 
times it approaches a Gaussian, as expected for a diffusion problem. From Eq. (|15J) we read 
oSpi(t) =G(l,t\k,0). 

If we choose a shock measure (|12|) as initial condition of a system which obeys Eq. (|13|) 
the form of the shock is preserved in time but due to the diffusion the state of the system 
will evolve into a superposition of shocks. Thus when performing an ensemble average the 
density profile is not a sharp step but smears out in time (see figure |2J as seen in Monte- 
Carlo (MC) simulations. Nevertheless a typical configuration of a single systems shows a 
sharp shock. 

Eq. (|13Jl directly gives the diffusion coefficient D s and the shock velocity v s : 

D s = 2Si + S 2 + S 3 + 2<5 4 

(16) 

v s = 5 3 + 25 A - 2Si - S 2 . 

The product measures to the left and to the right of the shock position are only possible 
for a special choice of boundary dynamics hi and &l-i- A possible choice is always those 
boundary dynamics which is the effect of reservoirs with densities pi respectively p 2 [26]: 




FIG. 1: Ensemble average of a diffusive shock measure (schematically): Due to superposition the 
density profile smears out although the form of the shock is preserved in each realization. 

One imagines that the lattice is extended by two additional sites at each boundary and one 
calculates the action of h-i and /i_ 2 {hh-i and hi) assuming that the sites are occupied 
according to the product measures in the bulk and determines the effective rates for the 
boundary action. Thus, n_i and sl x are substituted by p±, V-i and sZ± by (1 — p) - the 
other sites are handled accordingly. In this paper we always consider the limit L — ► oo such 
that the shock is assumed not to hit the boundaries. 

Given specific microscopic processes we test for the existence of diffusive shock measures 
by the following procedure: 

a) Check whether \p\) is a product measure of the system with periodic boundaries. 



b) Set up the boundary processes b\ such that 



^2 hm + bl ) \P*> P" 2 ' fc ) = c *-i IPi> P 2 ' fc ) 



(17) 



\m=l 



where Ck-i acts on sites k — 1 to k + 2. 

c) Check whether \p 2 ) is a product measure of the system with periodic boundaries. 

d) Set up the boundary processes &l-i such that 

^2 h m + 6 L _i J \p u p 2 , k) = 4+1 \Pl,P2, k) 

K m=k+1 / 

where dk+i acts on sites k + 1 to k + 2. 



(18) 



c) Check whether 



(c fc _i + 4+i) \pi,P2, k) = Sj. \pi,p2, k - 2) + 5 2 \pi, P2, k — 1) 

+ $3 \pi,p2, k + 1) + 84, \p!,p 2 , k + 2) 
- $s\Pi,P2,k) (19) 



In detail we do the following: After having set up the local Hamiltonians hk the action on 
a product measure is brought into a diagonal form as described above, 

hk \ P ) = C g \ P ) (20) 

where h k iag contains only the operators n m and 1 (v can be eliminated using v — 1 — n) 
acting on sites m = {k, k + 1, k + 2}. 

The requirement of a product measure (jsj) leads to the condition that the application 
of the Hamiltonian of the periodic system has to produce "telescope" -sums of diagonal 
operators which results in five equations of the rates for each density. 

Condition |gj) can be checked using the identity 

\p!,p 2 , k - 1) = ( -v k + —n k ) |pi, pa, k) (21) 

V - Pi Pi J 

and analogous identities for |p 1; p 2 , k — 2), |p 1; p 2) fe + 1) an d |/?i, P21 k + 2). A comparison of 
coefficients of the (only diagonal) operators on the left and right hand site of Eq. (|T9*j) then 
leads to another seven equations of the rates, pi (i = 1, 2) and Sj (j — 1, . . . , 5). Conservation 

of probability implies 

4 

*b = $>' (22) 

i=l 

which is not an additional equation but simplifies the calculation. 

If one of the densities is either zero or one, the Hamiltonian cannot be brought to diagonal 
form using Eqs. (jlUj) and (|21|). But in this case the action of the Hamiltonian simplifies and 
an analogous comparison of coefficients of creation and annihilation operators is possible. 



IV. CLASSIFICATION OF MODELS 

In order to find the models with three-site interactions which exhibit shock diffusion one 
could in principle try to solve Eq. (JT3j) for the general Hamiltonian. But as there are 56 
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microscopic processes (transitions from any of the 2 3 = 8 states to any different state) this 
task is tedious. 

We facilitate the procedure by writing a computer program which does all the symbolical 
calculation (transformation of the Hamiltonian into diagonal form, gathering of coefficients) 
and sets up the constituting set of equations. These equations are then solved by standard 
mathematical software. Still, the general solution of the problem is far too complex to extract 
useful informations. It is hence useful to investigate physically motivated sub-classes, as 
done in the following. 

A. Shocks from p\ = to p 2 = 1 

The completely empty lattice, p = 0, and the fully occupied one, p — 1, are nonfluctuating 
states. They are the two ground states of the zero-temperature Ising model. In this case only 
few processes play a role. All processes starting from 000 or AAA are forbidden because 
otherwise product measures with densities and 1 would not be stationary solutions of the 
system. The processes which act on the configurations 0A0, A00, A0A and AA0 are 
ineffectual because none of these configurations is possible, if the system is initialized with 
a shock measure. The remaining processes are 



00,4 ^ AAA; 00 A ^i 0AA; 00 A ^ 000; 
0AA ^ AAA: 0AA ^ 00 A: 0AA ^ 000. 



(23) 



The diffusion constants are then given by 

5 1 = K - 8 2 = K 1 + K 3 ; 5 3 = K 2 + K 4 ; 5 4 = K b . (24) 

The solution for a shock between p± = 1 and p 2 = is obtained by exchanging the roles of 
particles and holes (A «-► 0). 

This dynamics is a generalized zero temperature Ising model where the domain wall 
between spin up and spin down regions diffuses freely. As the space symmetric processes 
have no influence on the upward shock they can be included giving rise to a model allowing 
for both, upward and downward shocks. All ten rates are independent such that the drift 
of the domain walls can be chosen freely. 



B. Shocks from p x £ (0; 1) to p 2 G (0; 1) 

For each Hamiltonian with three-site interactions there is either none, exactly one product 
measure with a density in the open interval p G (0; 1), or infinitely many The latter is only 
possible for particle-conserving Hamiltonians. The proof of this assertion is possible for the 
case of three-site interactions, but rather technical and is therefore omitted here. It is based 
on analyzing the conditions on the rates and densities for the existence of a product measure 
Eq. (|11J1 for the general three-site interaction Hamiltonian. 

Instead we present a general consideration why the existence of two fluctuating station- 
ary product measures is impossible for non particle-conserving Hamiltonians. This is not 
intended to be a mathematical rigorous proof, it is rather a physical consistency check. This 
consideration is related to the positive rates conjecture, see Refs. J22J, 



23). 



Let us assume that there are two fluctuating stationary densities p\ ^ p 2 for a non 
particle-conserving Hamiltonian. Then due to fluctuations there is a finite probability that 
in the stationary state with density p x a region of density p 2 emerges. As p% is a stationary 
density this region has to be suppressed and vanishes again. Consequently the domain walls 
Pi\p2 and p 2 \p± are moving toward each other shortening the domain of p 2 . But if this is the 
case p 2 cannot be a stationary state, as a region of p\ emerging in a phase of p 2 is growing. 
Hence the assumption of the existence of two fluctuating stationary densities was wrong. 



Note that Gacs error correcting model |22t l23j | does not represent a counter example 
for this statement. It is crucial that product measures are considered: For this case no 
dynamics is able to determine the length of a one-dimensional domain - the domain growth 
is independent of the domain size, while phase separation requires a faster growth of larger 
domains. 

In Fig. El the time evolution of a non particle-conserving system is shown. The particle 
density in the low density region is not stationary and increases until the first shock vanishes. 
The second shock disappears due to another mechanism, here the diffusion equation of the 
shock is not fulfilled and consequently the shock dissolves. 

The situation is basically different if particle-conserving Hamiltonians are considered. In 
this case the system is not ergodic and hence a region of different density cannot evolve in 
the bulk of the system. Consequently several stationary densities are possible. 

As non-conserving Hamiltonians do not allow for the existence of two fluctuating phases 
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FIG. 2: Monte Carlo simulation of the model A00 — > AAA; AAA — ► 00 A, the number of sites 
is 500 and an average over 5000 systems was performed, the boundary conditions are periodic. As 
initial condition we took a shock measure with steps at site 166 from 0.01 to 1/2 and at site 333 
from 1/2 to 0. The left step is not stable because the product measure with p\ = 0.01 is not a 
stationary state and the right step is not stable because there the diffusion equation is not fulfilled. 



we only need to investigate particle-conserving Hamiltonians and recover the solution of 
the investigation of nearest neighbor interactions, the asymmetric simple exclusion process 
(ASEP) U: 



<*i = 0; 



&> 



A0 -5 0A; 
Ur- 

1-Pl 

1-P1 



--D, 



0A -> A0 

P2 



1 - 92 



(25) 



D L ; 



x l ~ P2 n 



S d = 



1 - p 2 " I- pi 

The inclusion of next-nearest neighbor interaction does not lead to further models in this 
case. 



C. Shocks from p\ = to p 2 € (0; 1) 

While the existence of two fluctuating stationary states is impossible, the existence of 
one fluctuating and one non fluctuating phase is easy to construct. The density p\ = is 
stationary if all processes from the empty lattice, 000, are prohibited. This is a violation of 
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the assumption underlying the positive rates conjecture as certain types of fluctuations are 
absent. Therefore two stationary states are possible, however, the nonfluctuating phase is 
unstable and the fluctuating phase will always enlarge at the expense of the nonfluctuating 
one. 

For the existence of a product measure p<i G (0; 1) it is then necessary that as well no 
process to the empty lattice is present. Apart from these constraints no further processes 
can be excluded a priori. We restrict ourselves to the case of a shock between p\ — to 
P2 *= (0; 1). The cases of a shock from p\ = 1 to p 2 G (0; 1) can be obtained by exchanging 
particles and holes, and the case of the fluctuating phase to the left can be obtained by a 
parity transformation. 

As argued above, the degrees of freedom are too many for a complete investigation and 
we classify the systems by the symmetries charge (C), parity (P) and time (T). In this 
context, charge symmetry is the invariance of the microscopic processes under the exchange 
of particles and holes (A <-► 0), i.e., for each process in the model there exists the C- 
symmetric one with the same rate. Although this picture is quite artificial when applied 
to particles it is natural in the language of spins where in the absence of an external field 
the symmetry between the up and down spin is obvious. Parity symmetry is the invariance 
of the microscopic processes under the exchange of left and right, i.e., for each process 
in the model there exists the P-symmetric one with the same rate. Time symmetry is the 
invariance under the transformation of t —* —t. A stochastic model is T-invariant if detailed 
balance with respect to its stationary distribution is fulfilled [25], i.e., these models are able 
to reach an equilibrium steady state. This is the case if for all states 771,772 the transition rates 
w m -* m , w V2 -> m and probabilities of finding the system in the configurations P(r)i), P(rj 2 ), 
obey the equation w m ^ m P{rji) = w V2 -> m P(r]2), i.e., there is no net current between states. 
If the stationary state of a system is a product measure, the validity of detailed balance is 
easily checked or refuted, because the calculation of the probabilities of the system to be in 
a specific states is trivial and additionally the absence of correlations permits to investigate 
only the local rates instead of the configurations of the whole lattice. 

One has to distinguish between a combined symmetry, for example PT, from an indepen- 
dent symmetry, here expressed by the symbol 'A', for example PAT. While in the former 
case the system is invariant after applying the symmetries one after each other, in the latter 
case the system is invariant under the symmetries applied each by themselves. 

12 



1. C A P AT symmetric systems with shock measures 

The 56 possible three-site interaction transitions can be arranged in 11 minimal models 
which obey C A P A T symmetry each (but are not necessarily described by exact diffusive 
shock measures). If we exclude those which have transitions involving the configurations 
000 or AAA only 6 remain. There are 63 combinations which can be build out of 6 
elements. We checked all combinations and found that there is only one model which is 
C A P A T symmetric and has an exact diffusive shock measure as solution: 

A0 i 0A; 0A -^ A0; 

00 A -^ 0AA; 0AA -4 00 A\ A00 -^ AA0; AA0 -4 A00; (26) 

,02 = 1/2; 5 1= 5 4 = 0; 5 2 = 2; 6 3 = 1. 

In this model diffusion is combined with branching processes which are only possible to 
a neighboring lattice site if the subsequent site is empty and its reversal, the coalescence 
process. 

Due to the CAP symmetry this model has the property that both the shock from or 1 
to 1/2 and the shock from 1/2 to or 1 is stable. Hence double shocks 0|||0 and 0|||1 are 
possible. We calculated the time evolution of an initial double shock (0|||1) in a periodic 
system by a Monte Carlo simulation as shown in Fig. El As calculated the shock moves with 
a drift of 1 and such that always the fluctuating phase penetrates the non-fluctuating phase. 
The shock 1|0 between site 500 and evolves into two steps, 1|| and ||0, since these shocks 
drift into the nonfluctuating phase and split up. As a consequence of the periodic boundary 
conditions the two shock fronts moving into the non-fluctuating phase coalesce after finite 
time and the non-fluctuating phase vanishes. This is a consequence of the stability of the 
fluctuating phase as argued in section IIVBI A detailed discussion of the reaction of shock 
fronts will be given inlVl 

2. CAP symmetric systems with shock measures 

There are 17 minimal models obeying CAP symmetry, if we exclude those which have 
transitions involving the configurations 000 and AAA only 8 remain. Thus in this case 
255 combinations have to be checked, but no additional model besides the C A P A T model 
can be found. 

13 
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FIG. 3: Monte Carlo simulation of the C A P A T symmetric model, the number of sites is 500 and 
an average over 10000 systems was performed. As initial condition we took a shock measure with 
steps at site 166 from to 1/2 and at site 333 from 1/2 to 1. 

3. P AT symmetric systems 

There are 18 minimal models obeying PAT symmetry, if we exclude those which have 
transitions involving the configurations 000 only 13 remain. Thus in this case 8191 com- 
binations have to be checked. We found 14 models, as presented in the appendix. Due to 
the P symmetry in each of these models downward shocks are also stable. 

In Fig. 0] we show a Monte Carlo simulation of the model D (see appendix) with u> — 1/2 
on a ring. As predicted two aspects can be observed: First, both the upward and the 
downward shock are stable. Second, the fluctuating phase spreads into the nonfluctuating 
one until the latter finally vanishes. We remind that due to the superposition of shocks the 
ensemble average does not exhibit a sharp step although each single realization does. 

4- C AT symmetric systems 

There are 16 minimal models obeying C AT symmetry, if we exclude those which include 
the configurations 000 or AAA only 9 remain. Thus in this case 511 combinations have 



14 



P 0.5 



i ' i ' r 



t=0 



i^ov»^».,» v 



J i L 



_L 



i ' r 



t=40 




l|W»|»fM|l^ll<#t»« 



J i L 



P 0.5 



100 200 300 400 500 100 200 300 400 500 

1 



"i ' r 



~T 



t=80 



r 



♦..sw.,.1. .,,.,!«< ..t.i^i.i,. IM rHK > 



% 



J i L 



"i ' r 



t=120 



V^i|..*4,'ltl»Wlf't'*">fl " »**W-' "H '«~> " »' " W W 



J i L 



0.5 



100 200 300 400 500 100 200 300 400 500 
k k 



FIG. 4: Monte Carlo simulation of the PAT symmetric model D with to = 1/2, the number of 
sites is 500 and an average over 5000 systems was performed. As initial condition we took a shock 
measure with steps at site 166 from to 1/3 and at site 333 from 1/3 to 0. 

to be checked. Besides the C A P A T model we find the following two models: 

0A0 -+ 0AA; 0AA -+ 0A0: A00 -^ A0A; A0A -+ A00: 

(27) 
p 2 = 1/2; 5 l = 5 2 = 5 3 = 5 i = 0. 

In this model only branching to the right is possible in the presence of another zero at the 
nearest neighboring site and the corresponding coalescence processes. The domain wall 0|| 
is not fluctuating even though the domain is. 

The second model is a combination of branching and coalescence processes to both di- 
rections 



00 A ->■ 0AA: 



A00 -> 0AA: 



0AA -+ 00 A; 

0AA -+ A00; 

P2 = 1/2; 



A00 ->• AA0; AA0 ->• A00; 
0A0 -^ A0A; A0A -+ 0A0; 
Si = 84 = 0; 5 2 = S 3 = 2. 



(2f 



The domain wall performs an unbiased diffusion with diffusion constant D = 2. 

For these two CAT models a downward shock has not necessarily to be stable as it is 
the case for the models which are P symmetric. Indeed, a downward shock is not stable in 
the two models because otherwise the space reflected versions of the processes would show 
a stable upward shock and would constitute additional CAT models. 
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5. Some further models 

Here we present some models which do not belong to the classes presented above. 
The following model is the totally asymmetric exclusion process combined with activated 
Langmuir kinetics used in Ref. [3J] to show hysteresis in driven diffusive systems: 

A0 -^ 0A; 
A0A ^ AAA; AAA % A0A; (29) 

p 2 = - — ; 5i = 5 2 = 5 4 = 0; 8 3 = (1 - p 2 ). 

The simplest model with a fluctuating shock front not included in nearest neighbor in- 
teraction models is 

A00 -► 00A; 5 4 = (1 - p 2 ) 2 , (30) 

where particles are only allowed to hop over a vacancy to the right. This model is PT 
symmetric, because reversing the direction of the process is equal to exchanging left and 
right. 

Another simple model is 



A00 ■!> AAA: AAA -^ 00 A\ 



L ! 



U 



_ :LZ — ^_ St = 62 = 8 3 = 0; 8 A = upj. 

uj — 1 



(31) 



This model is again PT symmetric by the same argument as in the model presented above. 
In this case one has to pay attention to the rates when reversing the time direction, the time 
reversed process is: 

AAA ^ A00; 00 A -i AAA. (32) 

For uj = 1 the density is p 2 = \- 

V. REACTIONS OF DOMAIN WALLS 

Up to now we have only investigated whether the shock fronts are stable and if so how 
they move. The movement of the shocks can be used to describe the dynamics of the systems 
as their positions are sufficient to characterize the state of the system. This has been used 

n 

for example in Ref. |3J where the dynamics with many degrees of freedom could be reduced 
to an effective one particle system. In order to describe a system completely by the position 
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of the domain walls we additionally have to investigate how the domain walls affect each 
other. 

A special case would be if the domain walls do not interact at all. This means that their 
rates do not change in the presence of another wall - certainly the possibility of mutual 
annihilation has to be included as a phase may vanish if two boundaries meet. In this case 
one could interpret the dynamics as the motion of annihilating random walkers. For those 
systems a direct link to free fermion system has been discussed in Ref. J4] and thus, for 
the unaffected movement of domain walls a description by free fermions could be possible. 
To this end one may apply the Jordan-Wigner [2j| transformation which converts spin- 1/2 
operators into fermionic creation and annihilation operators. 

The significance of the Jordan-Wigner transformation in this context is that some Hamil- 
tonians of spin-1/2 systems transform into fermion Hamiltonians which include only bilinear 
expressions of the fermionic operators - this can be regarded as a system of free fermions. 
For these Hamiltonians additional techniques for calculating the dynamical properties are 
available. Interestingly the dynamics of all particle models showing stable shock fronts with 
two-site interactions can be represented by the free motion of domain walls without in- 
teractions 14J. For the BCRW and the AKGP this is directly related to the free fermion 
character of these systems. 

It is the purpose of this section to investigate potential relations of the three-site interac- 
tions models found above to fermion systems, since a link between the domain wall motion 
and free fermion behavior would be interesting. We first discuss the dynamics of domain 
walls in detail and then turn to the transformation into fermion systems. 

6. Dynamics of domain walls 

When investigating the interaction of domain walls, the situation simplifies again by the 
fact that within the domains the probabilities are given by a product measure. By this two 
domain walls may only influence each other if the distance is smaller than three lattice sites 
and thus only operators acting on a small range have to be included. 

The first model to be considered is the C A P A T model. In Fig. El on the one hand it 
is shown how a domain wall 1|0 splits into two domain walls l\h and ||0 and on the other 
hand how two domain walls ||0 and 0|| coalesce. These two cases are now studied in detail 
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FIG. 5: The difference I — k of the two shock fronts in the state \k,l) perform a biased random 
walk. 

analytically. 

If the system is characterized simply by the position of the domain walls without further 
correlations it will be sufficient to describe the dynamics by states 



\k, Z; p x , p 2 , pz) = • • • ® \pi) ® • • • |pi) ® |p 2 ) ® • • • 

fe fc+i 

In order to describe the evolution of the step 1|0 we define 

\k,l) = \k,l; 1,^,0). 



IP2) ® |P3> ® ■ 
(-1 J 



(33) 



(34) 



In the following the densities are omitted for the sake of simplicity. By applying the Hamil- 
tonian of the C A P A T model on this state one gets 



|m = < 



4|A;-l,ZH-l)-4|fc,Z) l-k = l 

2\k-l,l) + 2\k,l + l) -4\k,l) l-k = 2 

2 |fc - 1, Z) + |ife + 1, Z) + 2 |fc, Z + 1) + \k, I - 1) - 6 \k, I) 

l-k>3. 



(35) 



Thus the time evolution of the system can be completely described by the movement of 
the domain walls, no additional correlations evolve. Starting from 1|0 both domain walls 
simultaneously move by one lattice site creating a fluctuating domain of two sites with 
density |. The time evolution of a state in which the domain walls are separated by a lattice 
site is given by the separate movement of both domain walls by one lattice site. A state 
in which the domain walls are separated by two or more lattice sites evolves simply by the 
rates 82 = 2 and £3 = 1 as calculated before. The time evolution of the distance of the 
domain walls Z — A; is illustrated in Fig. Once the domain walls are separated it is not 
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possible that they coalesce again as the / — k = 1 is an isolated point. Thus the movements 
of the domain walls are not independent, a repulsive interaction is present. 

Next, the interaction of the shock fronts ||0 and 0|~ shall be investigated. To this end 
we now define 

|fc,Z) = |M;i,0,i), (36) 

and test again whether the dynamics can be described in terms of the \k,l). Applying the 
Hamiltonian on the state \k — 1, k + 1) shows that additional correlations appear and that 
consequently the time evolution cannot be described simply by the location of the domain 
walls. 

Nevertheless it is instructive to analyze the appearance of the correlations in detail in 
order to reveal the link to free fermion systems. The correlations can be compensated by 
including the process A0A *=; AAA. By this the C symmetry is broken (the p — 1 phase is 
not stable anymore), but it is still a P AT model. 

Choosing the rate for the forward and backward reaction to be 2 one gets: 



!' M > 



l-k=l 

4 1|) + |ife - 1, 1) + \k, I + 1) - 6 \k, I) l-k = 2 

\k - 1,1) + \k,l + 1) + 2\k + 1,1) +2\k,l - 1) -6\k,l) 

l-k>3, 



(37) 



where ||) = \k, I) \k-i=i is the product measure with density p — |. This state is stationary 
which is recovered by the vanishing time derivative for k — I = 1. The time evolution of the 
system can be described completely by the movement of the shocks and their dynamics is 
independent of each other until they meet, then both are annihilated. 

Consequently this system is a candidate for the description by free fermions. However, 
by including the process A0A ^ AAA we get the model O of the appendix which is the 
BCRW - only a two-site interaction model. The transformation of the BCRW into a free 
fermion system is known (4J, [jjj . 

As a second example the model \j\ of the appendix is chosen as it is one of the simplest 
models. Again we define 

\k,l) = \k,l;p,0,p), (38) 
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and apply the Hamiltonian. The description by the states \k, I) is only closed if we set the 
parameter of the model w — 1 for which p = 1/2 and some of the rates vanish. One gets 



9 I/ ;\ 



o Z-fc=l 

4|i) + |fc-l,/) + |M + l)-6|M) l-k = 2 

\k - 1,1) + \k,l + 1) + 2\k + 1,1) + 2\k,l - 1) -6\k,l) 

l-k>3, 



(39) 



Next the time evolution of a state 



\k, l) = \k,l; 0,^,0) (40) 

is investigated. 
We find 

f l-k=l 

d 

qj. % = S |ife + 1, /) + |ib, I - 1) + 2 \k - 1, Z) + 2 |ife, Z + 1) - 6 |ife, Z) ( 41 ) 

l-k>3. 

For fc — Z = 2 it turns out that the time evolution of this state cannot be described by a 
superposition of shock measures, i.e. additional correlations emerge. Thus no independent 
movement of the shock fronts is possible in this case. This suggests that this is not a free 
fermion model. 



VI. CONCLUSIONS 

In this paper we have investigated exact diffusive shock measures in one-dimensional 
reaction diffusion systems with next nearest neighbor interactions and open boundaries. We 
distinguish the following three cases: 

1. The connection of two non-fluctuating phases, the two densities are and 1. The 
conditions that both states are stable exclude many models and we find many next nearest 
neighbor models as solution generalizing the Glauber Ising model at zero temperature. 

We restricted ourselves to the case of completely ordered connected phases as initial 
conditions. It would be interesting to investigate how a system evolves out of random initial 
conditions. In this scenario coarsening of the ordered domains or the emergence of a third 
stationary state which is fluctuating is possible. 
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2. The connection of two fluctuating phases, both densities are between and 1. It 
is argued that in general non-conservative models cannot have two fluctuating product 
measures as solution, in agreement with the positive rates conjecture. Consequently only 
conservative models have to be investigated and we recover the ASEP as the most general 
solution. Hence, the inclusion of three-site interactions does not lead to models not known 
from the investigation of nearest neighbor interactions. 

3. The connection of a non-fluctuating phase to a fluctuating phase, one density or 1 
and one between and 1. In this case numerous models exist and we classify the systems 
with respect to their symmetry. There is only one model which is C A P A T invariant, and 
this is as well the only model which is CAP invariant. Two additional models are found 
that are C AT invariant and 14 models are found that are PAT invariant. 

We stress that the mechanisms of exact shock measures are not suitable to construct a 
(one-species) model which shows phase separation on a ring. On the one hand, although 
models with one fluctuating phase and one non-fluctuating phase allow for stable up- and 
downward-shocks, the non-fluctuating phase will always vanish on a ring because the two 
shock fronts always enlarge the active region. On the other hand conservative models, which 
in principle allow for shocks between two fluctuating phases, are unable to show both up- and 
downward-shocks. This is a consequence of the collective velocity v c (p) which describes the 
movement of the center of mass of a disturbance in a region of a certain density p. In order 
that a shock is stable a disturbance has to tend toward the shock, v c (pi) > v s > v c {p2), where 
v s is the shock velocity. Obviously this equation can only hold either for the upward- or the 
downward-shock. Note that this consideration only holds for short-ranged, homogeneous 
one-species models, it is known that phase separation is possible in models with defects J28| 
or several species of particles J29|, |3(| . 

Although it is shown that double shocks 0|p|0 are possible we argue that this cannot be 
used to construct models which show phase separation in one-dimensional, short-ranged 
periodic systems with a single species. 

We have also investigated the influence of shock fronts on each other. A case of special 
interest is when the fronts move independently, except for the possibility of mutual anni- 
hilation. In the C A P A T model this cannot be observed, it turns out that it has to be 
combined with an additional process violating C symmetry. But by this the BCRW which 
is a two-site interaction model is recovered for which the independence of shock fronts is 
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known. 

We conclude that there is no direct connection of models whose time evolution is given 
by exact diffusive shock measures and free fermion systems. 

VII. APPENDIX: PAT MODELS 

The 14 P A T models are: 



Model A: 



A0A ^ AAA; 



AAA ^ A0A 



P 



;5i = 5 2 = S 3 = 8 4 







(42) 



00 a + 00 d 

The shock position is fixed without fluctuations in this model, but the model is 
not ergodic. 



Model B: 



UO 



A0 -+ 0A; 0A -+ A0\ 



00 A ^ 0AA; 0AA -+ 00 A; 

A00 ^ AA0; AA0 -+ A00; 

1 



(43) 



P2 



5 l = 5 4 = 0; 



So 



6. = 1 



00 + 1 " 1 - p 2 

In the case a; = 1 this model is C-invariant and we recover the C A P AT model. 



Model C: 



u> 



f>2 



00 + 1' 



A0 -^ 0A\ 0A -+ A0 



A0 ^ AA: 0A ^ AA: 



1 



5 1 = 5 4 = 0; 



1 



AA ->• A0; AA ->■ 0A; 

1 



(h 



I-P2 



(44) 



5, 



This model is a purely two-site interaction model and known as the branching 
coalescing random walk. It can be obtained by combining model A and B. 
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Model D: 



00 A ->• 0A0; 
0A0 -+ 00 A; 
00 A ™ 0AA; 
A00 ™ AA0; 
0A0 -^ A0A: 



A00 ->• 0A0; 
0A0 -+ A00; 
0AA A 00A; 
AA0 A A00; 
A0A -+ 0A0; 



0AA ^ AA0; AA0 ^ 0AA; 



0A0 ^ AAA; 
(l-2w) 



CJ 



P2 



W+l' 



1-W 

<5i = <5 4 = 0; 



1 
AAA ^ 0A0; 



cj < 



1 ~pi 



For w= i we get ct = and some of the rates vanish. 



Model E: 



Pi 



UJ 



u 



1' 



00A -^ 0A0; 
0A0 -+ 00A; 

00A ^ 0AA; 
A00 ^ AA0; 
0A0 ^ A0A; 
0AA -^ AA0; 
0AA ^ AAA; 
AAA -^ 0AA; 
5 1 = 5 4 = 0; 



A00 -^ 0A0; 
0A0 -^ A00; 



0AA -^ 00 A; 
AA0 -^ A00; 
A0A -^ 0A0; 
AA0 -^ 0AA; 
AA0 ^ AAA; 
AAA ±> AA0; 

s 2 - 2 • 



(45) 



(46) 



1 - Pi 
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Model F: 



00,4 ->• 040; 

040 -+ 004; 

040 aj(1 l > 2w) 0AA: 



0AA M" 040; 



u> 



Pi 



id 



1' 



0AA 2< ^i 1 AA0; 

9, ,2 
040 -^ AAA; 

1 
<Ji = S 4 = 0; 



For a; = h some rates vanish. 



A00 ->• 040; 
040 -+ 400; 
040 w(1 J^ 2a,) 440; 



440 



040; 



040 -^ 404; 404 -+ 040; 



440 ^ 044; 
444 -H- 040; 



1 -P2 



(47) 



Model G: 



P2 



CJ 



L)+ 1 



004 -+ 040; 
040 -+ 004; 
040 -A 044; 
044 -+ 004; 
040 A A0A; 
0AA -+ 440; 



044 ^ 444; 



444 A 044; 
5 1 = 5 4 = 0; 



400 -^ 040; 
040 -+ 400; 
040 ^ 440; 
440 A 040; 
404 -+ 040; 
440 A 044; 



440 ^ 444; 



444 A 440; 

2 



^ 



(48) 



1-/02 
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Model H: 



00,4 ->• 0A0; 
0A0 -+ 00 A; 

. 2ui(l— u>) . . 

0A0 U 0AA; 

. . 2(l-w) 

0AA — >■ 0A0: 



P2 



CO 



0AA ^ AA0; 
0A0 % AAA; 



w < l; 

<5i = 5 4 = 0; 



For u; = 1 some of the rates vanish 



A00 ->• 0A0; 
0A0 -+ A00; 
0A0 U AA0; 

. . 2(l-w) 

A40 — >■ 0A0; 



0AA -> A0A; A0A ->■ 0AA: 



A0A ->• AA0; 0AA ->• A0A; 



AA0 ^ 0AA: 



AAA ->• 0A0; 



1 ~P2 



(49) 



Model I: 



00 A ->• 0A0; 
0A0 -+ 00A; 
0A0 ^4 A0A: 



0A0 -> AAA; 
0AA v -» AAA; 



AAA 2 -^ 0AA; 



P2 



CO 



u < 2; 
^ = <5 4 = 0; 



CJ+ 1' 

For uj = 2 some of the rates vanish 



A00 ->• 0A0; 
0A0 -+ A00; 
A0A -+ 0A0; 



0AA ->■ AA0; AA0 -> 0AA; 



AAA A 0A0; 

w(2-w) 



AA0 



AAA; 



AAA 2 -4^ AA0; 



1 - P2 



(50) 
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Model J: 



LO 



P2 



A0 ->• 0A; 
0AA -^ AA0; 

0A0 -^ AAA; 



0AA ^ AAA: 



AAA ^ 0AA; 
00 > 1; 
<Si = £4 = 0; 



0A ->• A0; 
AA0 -^ 0AA; 

AAA A 0A0; 



AA0 "U 1 AAA: 



AAA ^ AA0: 



Model K: 



LU+ 1 

For ou — 1 some of the rates vanish. 

00A 
A00 
00A 
00A 
A0A 
00A 
A00 
00A 
AAA 
0A0 

1 r „ 

P2 = -; di = 4; 



1 ~P2 



^0AA; 
■^ AA0; 
-^ A00; 
-^ A0A; 
-^ 00A; 
-^ AA0; 
^0AA; 
-^ AAA; 
-^ 00A; 
-^ A0A; 

$2 = 2; 



0AA-^ 
AA0^ 
A00 -^ 
A00 -^ 
A0A^ 
0AA^> 
AA0-^ 
A00 -^ 
AAA^ 
A0A^> 

£3 = l; 



00A; 
A00; 
00A; 
A0A; 
A00; 
A00; 

0AA; 
AAA; 
A00; 
0A0; 

«5 4 = 1 



(51) 



(52) 
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Model L: 



Model M: 



P2 = ~; 



Pi 



00A 
00A 
A0A 
00A 
A00 
00A 
AAA 
0A0 
0AA 
0A0 
0AA 
AAA 
Si =4; 

00A 
00A 
A0A 
00A 
A00 
00A 
AAA 
0A0 
0AA 
0AA 
0AA 
0A0 

<*i=4; 



-+ A00; 
-^A0A; 
-+ 00 A\ 
^ AA0; 
-^0AA; 
^ AAA; 
-^ 00 A; 
^0AA; 
-^ 0A0; 
-^A0A; 
-^ AAA; 
-^0AA; 

$2 = 2; 

-+ A00; 
^ A0A; 
-+ 00i; 
^ AA0; 
^>0AA; 
-^ AAA; 
■1+ 00 A\ 
-^0AA; 
-^ 0A0; 
^0AA; 
-+ 00A; 

5 2 = 2; 



A00 -^ 
A00 -^ 
A0A^> 
0AA^ 
AA0-^ 
A00 -^ 

0A0 -^ 
AA0^> 
A0A^ 
AA0^ 

£3 = i; 

A00 -^ 
A00 -^ 
A0A^ 
0AA^ 
AA0^> 
A00 -^ 

0A0 -^ 
AA0^ 
0A0 -^ 

AAA-^ 
S 3 = l; 



00A; 
A0A; 
A00; 
A00; 
0AA; 
AAA; 
A00; 
AA0; 
0A0; 
0A0; 
AAA; 
AA0; 
5 A = 1 

00A; 
A0A; 
A00; 
A00; 

0AA; 
AAA; 
A00; 
AA0; 
0A0; 
AA0; 
0A0; 
0A0; 

5 4 = 1 



(53) 



(54) 
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Model N: 



P2 = ~; 



00A 
00A 
A0A 
00A 
A00 
00A 
AAA 
0A0 
0AA 
0A0 
Si =4; 



-^ A00; 
-^A0A; 
-^ 00i; 
^ AA0; 
^0AA; 
-^ AAA; 
■1+ 00 A\ 
-^ A0A; 
-^AA0; 



-^ AAA; 

$2 = 2; 



A00 -^ 
A00 -^ 
A0A^> 
0AA^> 
AA0^> 
A00 -^ 

AAA^ 
A0A-^ 
AA0-^ 
AAA-^ 

S3 = 1; 



00A; 
A0A; 
A00; 
A00; 
0AA; 
AAA; 
A00; 
0A0; 
0AA; 
0A0; 
5 A = 
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